Determination of elastomer material properties for the Mullins effect using a bi-axial test device

ABSTRACT

Systems and methods for determining material properties of elastomers for the Mullins effect are described. In one aspect of the present invention, material properties of an elastomer membrane specimen are obtained using a system comprising a bi-axial test device, a pump, a fluid reservoir, a linear variable differential transformer, a pressure transducer and a computer. The bi-axial test device comprises a top plate and a bottom plate. The top plate has a circular hole configured to allow the specimen to be expanded up by pressures of the inflating fluids. The circular hole is so dimensioned that the specimen can be expanded with a substantially low pressure. The bottom plate is a solid plate configured with a fluid intake at one side and a fluid outlet at the other end. The fluid intake is connected to the fluid reservoir. Fluids stored in the fluid reservoir are pumped into the bi-axial test device by the pump. The fluid outlet is connected to the pressure transducer.

FIELD OF THE INVENTION

The present invention generally relates to determination of rubber-like material properties, more particularly to methods and systems for determining the material properties of elastomers for the Mullins effect using a bi-axial test device.

BACKGROUND OF THE INVENTION

Rubber-like materials such as elastomers have been used in many parts and structures in various industries (e.g., automotive, aerospace, etc.) for years. But the mechanical properties (e.g., stress-strain or stress-stretch ratio relationship) of elastomer, other than the elastic properties, are still not clearly defined; hence, designs and analyses (especially used in computer-aided engineering) of these structures are generally based on the elastic properties of the elastomer only. In reality, elastomers exhibit non-elastic effects such as the Mullins effect, viscoelastic, and chronorheological behavior, and the magnitudes of the non-elastic properties are often large enough that they should not be neglected.

Elastomer in its virgin state exhibits a relatively stiffer response on the initial loading. When the elastomer is loaded, subsequently unloaded, then reloaded, the stress-strain relationship follows a significantly softer path. After several unloading-reloading cycles, the stress-strain relationship stabilizes, and additional unloading-reloading cycles retrace the stabilized path in the stress-strain curve. The non-elastic material behavior of elastomer described herein is referred to as the Mullins effect, in which the stress-strain relationship depends on the maximum loading previously encountered.

To date, very few analytical and experimental studies for determining non-elastic properties of elastomers have been attempted. This is because the study of mechanics for elastomers must consider both geometric and material nonlinearities. The additional effects and the lack of an adequate constitutive equation that describes these phenomena make the analytical and experimental studies very difficult.

The determination of the material properties is generally conducted using uni-axial test, as shown in FIG. 1A, in which an elastomer specimen 102 is pulled by a uni-axial tension 104 (i.e., one-dimensional test) at either end. However, the material properties obtained from the one-dimensional test may not represent true behaviors of elastomer contained in a structure, which is generally not in a one-dimensional space. To solve this deficiency, one of the prior art attempts is to stretch a sheet of an elastomer specimen (membrane) 112 in two planar directions as shown in FIG. 1B. The specimen 112 is pulled by equal tensions 114 in both of the specimens's planar axes. However, this prior art approach requires a huge piece of specimen (thereby huge laboratory) to avoid the edge effects from point loads (i.e., tensions 114) applied around the perimeter. This is not a practical solution because the test would have to be conducted in a very large test facility or laboratory.

Further, the current or prior art numerical equations (i.e., constitutive equations) are not adequate to represent true behaviors of elastomers. One of the prior art constitutive equations is the Ogden equation, which suggests or assumes the loading and subsequent reloading paths are the same. The Ogden equation does not represent the true behaviors of elastomers, which become softer in a subsequent reloading path than in the original loading path.

The elastomer behaviors as shown in FIG. 2 are calculated using the Ogden equation. In FIG. 2, the vertical axis represents stress in the elastomer. The horizontal axis represents the stretch ratio λ of the elastomer. The stretch ratio is defined as the stretched length divided by the original length of the elastomer. Therefore, the relationship between the stretch ratio λ and the strain ε is that ε=λ−1. The elastomer is first loaded following the path 202 starting at un-stretched position 201 (i.e., stretch ratio equal to one) until the stretch ratio reaches four (4) at 203. Then the elastomer is unloaded following a first unload path 204 back to the origin 201. Next the elastomer is reloaded following the reloading path 212 to a stretch ratio about 5.5 at 213. The first portion of the reloading path 212 is assumed to be the same as the original loading path 202 in the Ogden equation. It is emphasized that this is an incorrect assumption thus leading to inaccurate numerical analysis of elastomers. Finally the elastomer is unloaded again at 213. Based on the Ogden equation, the elastomer follows the second unloading path 214 back to the origin 201.

Given the foregoing drawbacks, problems and limitations of the prior art, it would be desirable to have improved methods and systems to determine material properties of elastomers such that the more correct behaviors of elastomers such as the Mulling effect can be numerically calculated and analyzed.

BRIEF SUMMARY OF THE INVENTION

This section is for the purpose of summarizing some aspects of the present invention and to briefly introduce some preferred embodiments. Simplifications or omissions in this section as well as in the abstract and the title herein may be made to avoid obscuring the purpose of the section. Such simplifications or omissions are not intended to limit the scope of the present invention.

The present invention discloses systems and methods to determine material properties of elastomers for the Mullins effect using a bi-axial test device. According to one aspect of the present invention, material properties of an elastomer membrane specimen are obtained using a material properties determination system including a bi-axial test device, a pump, a fluid reservoir, a linear variable differential transformer (LVDT), a pressure transducer and a computer. The bi-axial test device comprises a top plate and a bottom plate. The top plate has a circular hole configured to allow the specimen to be expanded up by pressures of the inflating fluids. The circular hole or opening is so dimensioned that the specimen can be expanded with a substantially low pressure (e.g., less than 15 pound-force per square inch (psi)). The bottom plate is a solid plate configured with a fluid intake at one side and a fluid outlet at the other end. The fluid intake is connected to the fluid reservoir. Fluids stored in the fluid reservoir are pumped into the bi-axial test device by the pump. The fluid outlet is connected to the pressure transducer.

The elastomer membrane specimen has a uniform thickness and a large enough size to be sandwiched between the top and bottom plates for forming a fluid-tight seal, when the top and bottom plates are coupled together with a plurality of connector means (e.g., screws). The amount of the pressure in the bi-axial test device is measured by the pressure transducer, while the inflated vertical displacement at the center of the specimen (i.e., center of the circular hole) is measured by the LVDT. Both the measured vertical displacement and the pressure are gathered and plotted in a computer. Since the pressure is uniformly distributed inside the inflated membrane specimen, the tensions at the center of the membrane specimen are equal in two planar directions (i.e., bi-axial tension). The size of the circular hole is relatively small with a radius between 1-2 inches.

In another aspect of the present invention, a new constitutive equation with separate loading, unloading and subsequent reloading paths is created. Two sets of material constants of the new constitutive equation for the unloading and reloading paths are calculated using least square fitting techniques based on the respective measured data from a bi-axial test. The new constitutive equation can be programmed in a computer aided engineering analysis software application such as finite element analysis software. As a result, the Mullins effect can be numerically analyzed for structure containing elastomers using the new constitutive equation including the measured data from the bi-axial test.

According to one embodiment, the present invention is a bi-axial test device for determining elastomer material properties for the Mullins effect under bi-axial tension, the device comprises at least the following: a top plate having a circular opening, the circular opening having a diameter; a bottom plate having an inflating fluid intake and an inflating fluid outlet; and a plurality of connecting means for coupling the top and bottom plate.

According to another embodiment, the present invention is a system for determining elastomer material properties for the Mullins effect under bi-axial tension, the system comprises at least the following: a bi-axial test device having an the inflating fluid intake and an inflating fluid outlet; a fluid reservoir coupling to the inflating fluid intake providing inflating fluids; and a pump, coupling to the fluid reservoir, configured to control the inflating fluids. The system further comprises the following: a linear variable differential transformer (LVDT) configured to measure vertical displacement at the center of an elastomer specimen placed in the bi-axial test device; a pressure transducer, coupling to the inflating fluid outlet, configured to measure pressure of the inflating fluids; and a computer, coupling to the LVDT and the pressure transducer, configured to gather the measured pressure and displacement.

According to yet another embodiment, the present invention is a method for numerically analyzing structure containing elastomer, the method comprises at least the following: conducting a bi-axial tension test of a specimen of the elastomer using a bi-axial test device of an elastomer material properties determination system, wherein the bi-axial tension test comprises loading, unloading and reloading phases; plotting measured pressure versus vertical displacement curve include the loading, unloading and reloading phases; determining a first set of material constants by a least square fitting technique using the measured pressures and displacements in the unloading phase; and determining a second set of material constants by least square fitting technique using the measured pressures and displacements in the reloading phase; and numerically analyzing the structure using a constitutive equation that uses the first and the second sets of material constants to describe behaviors of the elastomer.

One of the objects, features, and advantages of the present invention is to allow determination of elastomer material properties using the bi-axial test device in a relatively compact laboratory. Other objects, features, and advantages of the present invention will become apparent upon examining the following detailed description of an embodiment thereof, taken in conjunction with the attached drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

These and other features, aspects, and advantages of the present invention will be better understood with regard to the following description, appended claims, and accompanying drawings as follows:

FIG. 1A is a diagram showing an elastomer stick specimen under uni-axial tension;

FIG. 1B is a diagram showing a sheet of elastomer membrane specimen under bi-axial tension based on a prior art solution;

FIG. 2 is a diagram showing the stress-stretch ratio relationship of an elastomer under load-unload-reload phases calculated using the Ogden constitutive equation—a prior art approach;

FIG. 3A is a diagram depicting an exemplary elastomer material properties determination system, according to an embodiment of the present invention;

FIG. 3B is a diagram showing a lateral elevation view of the bi-axial test device of FIG. 3A and a top plan view of the elastomer membrane specimen in accordance with one embodiment of the present invention;

FIG. 3C is a diagram showing components of an exemplary bi-axial test device of FIG. 3A in accordance with one embodiment of the present invention;

FIG. 4A is a diagram showing an exemplary pressure versus vertical displacement curves generated by the elastomer material properties determination system of FIG. 3A;

FIG. 4B is a diagram showing the comparison between exact solution, approximate solution and test data of inflating a circular elastomer membrane based on Mooney material model.

FIG. 4C is a diagram illustrating a stress-stretch ratio relationship (curves) of an elastomer under load-unload-reload phases calculated using a new constitutive equation in accordance with one embodiment of the present invention; and

FIG. 5 is a flow chart showing the process of numerically analyzing structure containing elastomer with the material properties obtained from a bi-axial test using the bi-axial test device of FIG. 3A, according to an embodiment of the present invention.

DETAILED DESCRIPTION

In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will become obvious to those skilled in the art that the present invention may be practiced without these specific details. The descriptions and representations herein are the common means used by those experienced or skilled in the art to most effectively convey the substance of their work to others skilled in the art. In other instances, well-known methods, procedures, components, and circuitry have not been described in detail to avoid unnecessarily obscuring aspects of the present invention.

Reference herein to “one embodiment” or “an embodiment” means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the invention. The appearances of the phrase “in one embodiment” in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. Further, the order of blocks in process flowcharts or diagrams representing one or more embodiments of the invention do not inherently indicate any particular order nor imply any limitations in the invention.

Embodiments of the present invention are discussed herein with reference to FIGS. 3A-5. However, those skilled in the art will readily appreciate that the detailed description given herein with respect to these figure is for explanatory purposes as the invention extends beyond these limited embodiments.

Referring now to the drawings, in which like numerals refer to like parts throughout several views. FIG. 3A depicts an exemplary elastomer material properties determination system 30, according to an embodiment of the present invention. The system 30 comprises an inflating fluid subsystem and a data measurement subsystem, both coupling to a bi-axial test device 310. The bi-axial test device 310 comprises a top plate 314 and bottom plate 316. The top plate 314 includes a circular hole. The bottom plate 316 is a solid plate coupled with an inflating fluid intake 312 and with an inflating fluid outlet 318. An elastomer membrane specimen 315 is placed between the top plate 314 and the bottom plate 316 during a bi-axial tension test. The inflating fluid subsystem comprises a fluid reservoir 304 and a pump 302. The pump 302 is used for pumping inflating fluids (e.g., air or water), from the fluid reservoir 304 to the bi-axial test device 310, to inflate, deflate and re-inflate the elastomer membrane specimen 315. The data measurement subsystem includes a computer 326, a linear variable differential transformer (LVDT) 324 and a pressure transducer 322. The pressure transducer 322, coupling to the bi-axial test device 310 via the inflating fluid outlet 318, is configured to measure pressures of the inflating fluids throughout a bi-axial test. The LVDT 324 is configured to measure vertical displacement at the center of the elastomer membrane specimen 315 during the bi-axial test. Both the pressure transducer 322 and the LVDT 324 are coupled to the computer 326 (e.g., personal computer, server, laptop, desktop), which gathers and plots the measured data (i.e., pressure and displacement). The pump 302 is optionally configured to be controlled by the computer 326.

FIG. 3B shows a lateral elevation view 31 of the bi-axial test device 310 and a top view 32 of the elastomer membrane specimen 315 of FIG. 3A, according to an embodiment of the present invention. The elastomer specimen 315 has a uniform thickness H 332. As shown in the lateral view 31, the elastomer membrane specimen 315 is expanded up to a vertical displacement Δ 336 at the center 342 by a pressure 338 from the inflating fluids. The circular hole or opening of the top plate 314 has a radius R 334. The top view 32 shows that the membrane 315 at the center 342 is subjected to equal bi-axial tensions 344 as the elastomer membrane 315 is being expanded upwards by the pressure of the inflating fluids.

Sometime, a stiffer material will require higher inflating pressure, which may present a technical problem to achieve the higher pressure in a laboratory. The simplest way to solve this problem is to use either a larger size or a thinner specimen, for example, a membrane, doubled the size of radius or halved the thickness, will require one-half of the pressure to achieve a same vertical displacement. In one embodiment, typical size of the circular hole is relatively small with a radius between one to two inches.

The plane view of the top plate 314, the bottom plate 316 and the elastomer membrane specimen 315 (virgin state) are shown in FIG. 3C, according to an embodiment of the present invention. To conduct a bi-axial test, the elastomer membrane specimen 315 is sandwiched between the top plate 314 and the bottom plate 316. The top and bottom plates are securely coupled together with a plurality of screws. Each of the screws is screwed in a pair of screw holes 354, one in the top plate, the other the bottom plate. The circular hole or opening is so dimensioned that the specimen 315 can be expanded with a substantially low pressure (e.g., maximum of 15 psi).

In order to create a fluid-tight requirement for the bi-axial test, an O-ring is used for sealing the perimeter of the circular hole or opening. In one embodiment, two O-rings 352 are located on the top side of the bottom plate 316, while one corresponding O-ring (not shown) is located on the bottom side of the top plate 314. The elastomer membrane specimen 315 is so dimensioned that the size is large enough to cover the circular hole of the top plate 314, and to be able to form a fluid-tight environment with the top and the bottom plates. The elastomer membrane specimen 315 has a uniform thickness such that the inflating fluids inflate the specimen 315 uniformly. It is noted that shape of the top and bottom plates, number of screws, type of screws and O-rings are not limited to the embodiment as shown in FIG. 3C. Other means that can achieve the same purpose of providing a fluid-tight test environment with a circular hole may be used for the present invention.

Further, the material properties of elastomers are temperature-dependent in general. Therefore the Mullins effect is also temperature-dependent. The bi-axial test device 310, according to one embodiment of the present invention, is small enough such that the entire device 310 can be placed in a water bath or in an environment control chamber. The prior art approaches need huge specimen hence requiring huge temperature control facility.

FIG. 4A shows an exemplary pressure versus vertical displacement (P-Δ) curves 40 generated by the elastomer material properties determination system 30 of FIG. 3A. The vertical axis of the P-Δ curves 40 represents the pressure P 338, while the horizontal axis represents the vertical displacement Δ 336.

Using the elastomer material properties determination system 30, an elastomer membrane specimen 315 is inflated and measured. In the beginning of the bi-axial test, there is neither pressure nor displacement at 401. As more fluids are pumped into the bi-axial test device 310, the displacement and the pressure increase accordingly following path 402 to reach a first predetermined displacement (e.g., 0.8 inches) at 403. Then the specimen 315 is unloaded by reducing the pressure of the inflating fluids. The unloading phase following path 404 back to the original un-inflated state at 401. The specimen 315 is reloaded back to 403. The bi-axial test results show clearly that the reloading phase follows another path 412 a, which is softer than the original loading phase path 402 but stiffer than the unloading path 404. Two more cycles of unloading and reloading of the specimen 315 are conducted thereafter.

Next, at 403 after three cycles of loading and reloading, the inflating fluids are increased to inflate the specimen 315 to a larger or second pre-determined vertical displacement (e.g., 1.6 inches) at 413 following path 412 b. The path 412 b is a path of loading specimen 315 in its virgin state again. The bi-axial test continues after that by repetitively unloading the specimen 315 to the original state 401 and reloading back to the second vertical displacement at 413. It is clear that the unloading phase follows path 414 and reloading phase follows path 422. The P-Δ curves 40 clearly demonstrate the Mullins effect, which includes a softer reloading path after the initial loading and reloading.

In order to use the measured data from the bi-axial test in a numerical analysis, a new constitutive equation is created to describe the elastomer behaviors for the Mullins effect demonstrated by the P-Δ curves 40 in FIG. 4A. The new constitutive equation (i.e., strain-energy density function) including Mullins effect damage function of an elastomer is defined as {tilde under (W)}(λ_(i)). {tilde under (W)}(λ_(i))=ηW(λ_(i))  (1) where W(λ_(i)) is the strain-energy density function based on the initial loading (i.e., virgin state), η is a damage function for the Mullins effect, and λ_(i) is the stretch ratio of the elastomer.

For initial loading η=1  (2a)

For unloading

$\begin{matrix} {\eta = {1 - {\frac{1}{r_{1}}{\tanh\left\lbrack {\frac{1}{m_{1}}\left( {W_{m} - W} \right)} \right\rbrack}}}} & \left( {2b} \right) \end{matrix}$

For subsequent reloading

$\begin{matrix} {\eta = {1 - {\frac{1}{r_{2}}{\tanh\left\lbrack {\frac{1}{m_{2}}\left( {W_{m} - W} \right)} \right\rbrack}}}} & \left( {2c} \right) \end{matrix}$ W_(m)(λ_(i)) is the previous maximum of the strain-energy density function before unloading (e.g., 403 and 413 of FIG. 4A). r₁, r₂, m₁ and m₂ are the material constants for Mullins effect damage function. With the Mullins effect damage function, the loading 402, 412 b, unloading 404, 414 and subsequent reloading 412 a, 422 follow different paths as shown in FIG. 4A. It is noted that the material constants, r₁ and r₂, are greater than one in general.

Generally, the strain-energy function W(λ_(i)) may be in one of the many material model forms, for example, neo-Hookean, Mooney, Ogden incompressible or Ogden compressible. When Ogden compressible material model is considered, all other are special cases. For highly compressible materials the strain-energy density for Ogden constitutive equation can be written as follows:

$\begin{matrix} {W = {\underset{j = 1}{\sum\limits^{m}}{\frac{C_{j}}{b_{j}}\left\lbrack {\lambda_{1}^{b_{j}} + \lambda_{2}^{b_{j}} + \lambda_{3}^{b_{j}} - 3 + {\frac{1}{n}\left( {J^{- {nb}_{j}} - 1} \right)}} \right\rbrack}}} & (3) \end{matrix}$ where C_(j), b_(j) and n are material constants, J=λ₁λ₂λ₃ represents the ratio of deformed to un-deformed volume, and λ₁, λ₂ and λ₃ are stretch ratios of the elastomer in three orthogonal directions, respectively. When n=∞, Equation (3) becomes the Ogden incompressible material model. When n=∞, m=2, b₁=2 and b₂=−2, Equation (3) becomes the Mooney material model. When n=∞, m=1 and b₁=2, Equation (3) becomes the neo-Hookean material model.

According to one embodiment of the present invention, the neo-Hookean material model is used to demonstrate how material constants, r₁, r₂, m₁ and m₂, are determined from the measured data of the bi-axial test. The strain-energy function for the neo-Hookean material model is listed as follows:

$\begin{matrix} {{W\left( \lambda_{i} \right)} = {C_{1}\left( {\lambda_{1}^{2} + \lambda_{2}^{2} + \lambda_{3}^{2} - 3} \right)}} & (4) \end{matrix}$ An approximate relationship between the inflating pressure P 338 and the deformation at the pole Δ 336 can be obtained as follows:

$\begin{matrix} {P = {\frac{2C_{1}H}{R}\left\{ {\frac{1}{\left( {\frac{\Delta}{R} + \frac{R}{\Delta}} \right)}\left\lbrack {1 - \frac{1}{\lambda^{6}}} \right\rbrack} \right\}}} & (5) \end{matrix}$ The relationship between λ and Δ is

$\begin{matrix} {\lambda = {\frac{\left( {\frac{\Delta}{R} + \frac{R}{\Delta}} \right)}{2}{\sin^{- 1}\left( \frac{2}{\left( {\frac{\Delta}{R} + \frac{R}{\Delta}} \right)} \right)}}} & (6) \end{matrix}$ Where λ is the stretch ratio of the elastomer in the planar direction, R 334 is the radius of the circular membrane (e.g., elastomer specimen 315) and H 332 is the thickness of the circular membrane. Similar expression may be obtained for Mooney material model and other constitutive equations. The exact numerical solution 432, approximate solution 434 obtained from Equations (5) and (6), and test data 436 for inflating a circular thin disk (e.g., elastomer specimen 315) are compared and plotted in FIG. 4B. It is shown that the approximate solution 434 from Equations (5) and (6) represents a very good correlation to the true solution 432 and measured test data 436. It is noted that the material used in FIG. 4B is for Mooney material model with α=C₂/C₁=0.02.

When the Mullins effect damage function is considered, Equation (5) becomes

$\begin{matrix} {P = {\frac{2C_{1}H}{R}\left\{ {\frac{1}{\left( {\frac{\Delta}{R} + \frac{R}{\Delta}} \right)}\left\lbrack {1 - \frac{1}{\lambda^{6}}} \right\rbrack} \right\}*\eta}} & (7) \end{matrix}$ Using the measured data from the bi-axial test and Equation (7), the material constants in Equations (2b) and (2c) can be obtained. The first set of material constants, r₁ and m₁, are obtained from the pressure-deformation unloading curve (e.g., path 404 or 414 of FIG. 4A). The second set of material constants, r₂ and m₂, are obtained from the pressure-deformation reloading curve (path 412 a or 422 of FIG. 4A). These constants are all obtained from the least-square curve fit. Equations 2(a), 2(b) and 2(c) can be applied to one-, two- and three-dimensional analyses of structures containing elastomers tested.

Referring to FIG. 4C, which shows stress-stretch ratio (σ-λ) curves 42 based on Equations (2a), (2b) and (2c) in accordance with one embodiment of the present invention. The σ-λ curves 42 show similar behaviors with the measured test data (i.e., P-Δ curves 40) as shown in FIG. 4A. The initial stress σ and the initial stretch ratio λ (i.e., un-stretched state) start at the origin 441. A loading path 442 leads to a state 443. The stress σ follows a first unloading path 444 back to the origin 441. Reloading phase follows paths 452 a and 452 b to reach another state 453. Then the stress σ follows a second unloading path 454 to the origin 441. The subsequent reloading phase follows a path 462. The similarity of the curves 42 and the P-Δ curves 40 validates the usage of Equations (2a), (2b) and (2c) to numerically analyze elastomers.

FIG. 5 is a flowchart illustrating a process 50 for numerically analyzing structure containing elastomer with the material properties obtained from a bi-axial test using the bi-axial test device of FIG. 3A, according to an embodiment of the present invention. The process 50 is preferably understood in conjunction with the previous figures especially FIG. 3A and FIG. 4A, and the process 50 is preferably implemented in software.

The process 50 starts, at 502, by conducting a bi-axial test of an elastomer membrane specimen 315 of interest using an elastomer material properties determination system (e.g., system 30). Pressure versus vertical displacement data are measured and gathered following a load-unload-reload test (i.e., bi-axial test) described in FIG. 4A at 504. Using the geometry of the elastomer membrane specimen 315 (i.e., radius R 334 and thickness H 332), a first set of the material constants, r₁ and m₁, of Equation 2(b) are determined at 506. The determination is based on the measured data in the unloading phase (e.g., path 404 or 414 of FIG. 4A) using a least square fit technique. Similarly, a second set of material constants, r₂ and m₂, of Equation 2(c) are determined at 508. The difference is that the material data in the reloading phase (e.g., path 412 a or 422 of FIG. 4A) are used. After both the first set and the second set of material constants are determined, at 510, any structure containing the elastomer of interest (i.e., the tested elastomer at 502) can be numerically analyzed using constitutive equations (e.g., Equations 2(a), 2(b) and 2(c)). One implementation of this technique is to include the constitutive equation in a finite element analysis program, such that the user can model the structure with finite element with material properties defined by Equations 2(a)-2(c) with the first and second set of material constants defined automatically by the application program or manually by users.

The methods or algorithms described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executable by a processor, or in a combination of both, in the form of control logic, programming instructions, or other directions, and may be contained in a single device or distributed across multiple devices. A software module may reside in Random Access Memory (RAM), flash memory, Read-Only Memory (ROM), Erasable Programmable Read-Only Memory (EPROM), Electronically Erasable Programmable Read-Only Memory (EEPROM), registers, hard disk, a removable disk, Compact Disc Read-Only Memory (CD-ROM), Digital Video Disc (DVD) or any other form of storage medium known in the art. A storage medium may be coupled to the processor such that the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium may be integral to the processor.

Although the present invention has been described with reference to specific embodiments thereof, these embodiments are merely illustrative, and not restrictive of, the present invention. Various modifications or changes to the specifically disclosed exemplary embodiments will be suggested to persons skilled in the art. For example, whereas the bi-axial test device has been shown and described to determine the material properties of elastomers, other types of rubber-like material may be determined. Whereas the LVDT has been described and shown to measure the vertical displacement of the elastomer membrane, other types of measurement device ay also be used instead. Whereas the pressure transducer has been shown and described to measure the pressure in the bi-axial test device, other types of pressure gauge may be used. In summary, the scope of the invention should not be restricted to the specific exemplary embodiments disclosed herein, and all modifications that are readily suggested to those of ordinary skill in the art should be included within the spirit and purview of this application and scope of the appended claims. 

1. A device for determining elastomer material properties including the Mullins effect under bi-axial tension comprising: a top plate having a circular opening, the circular opening having a diameter; a bottom plate having an inflating fluid intake and an inflating fluid outlet; an elastomer specimen configured to be sandwiched between the top and bottom plate; and a plurality of connector means for coupling the top and bottom plates; wherein the elastomer specimen's area is larger than the circular opening so that a fluid-tight environment is formed when the top plate and the bottom plate are connected together by the connector means and wherein the fluid-tight environment further includes one first O-ring fixedly attached to bottom of the top plate and two second O-rings fixedly attached to top of the bottom plate, the first O-ring fits in between the two second O-rings.
 2. The device of claim 1, wherein the elastomer specimen is an elastomer membrane with uniform thickness.
 3. The device of claim 2, wherein the elastomer specimen is configured to be expanded upwards by pressures of inflating fluids.
 4. The device of claim 3, wherein the center of the elastomer membrane is under equal planar tension forces when the inflating fluids are pumped into the device.
 5. The device of claim 1, wherein the connector means are screws.
 6. A device for determining elastomer material properties including the Mullins effect under bi-axial tension comprising: a top plate having a circular opening, the circular opening having a diameter; a bottom plate having an inflating fluid intake and an inflating fluid outlet; an elastomer specimen configured to be sandwiched between the top and bottom plate; and a plurality of connector means for coupling the top and bottom plate; wherein the elastomer specimen's area is larger than the circular opening so that a fluid-tight environment is formed when the top plate and the bottom plate are connected together by the connector means and wherein the fluid-tight environment further includes two first O-rings fixedly attached to bottom of the top plate and one second O-ring fixedly attached to top of the bottom plate, the second O-ring fits in between the two first O-rings.
 7. The device of claim 6, wherein the elastomer specimen is an elastomer membrane with uniform thickness.
 8. The device of claim 7, wherein the elastomer specimen is configured to be expanded upwards by pressures of inflating fluids.
 9. The device of claim 8, wherein the center of the elastomer membrane is under equal planar tension forces when the inflating fluids are pumped into the device.
 10. The device of claim 6, wherein the connector means are screws.
 11. A system for determining elastomer material properties including the Mullins effect under bi-axial tension comprising: a bi-axial test device having an inflating fluid intake and an inflating fluid outlet; a top plate having a round opening, the round opening having a substantially circular shape; a bottom plate, wherein the inflating fluid intake and the inflating fluid outlet are operably coupling to the bottom plate; and a plurality of connector means for coupling the top and bottom plate; wherein an elastomer specimen is configured to be sandwiched between the top and bottom plate; a fluid reservoir, coupling to the inflating fluid intake, configured to store inflating fluids; a pump, coupling to the fluid reservoir, configured to control the inflating fluids; a linear variable differential transformer (LVDT) configured to measure vertical displacement at the center of the elastomer membrane specimen placed in the bi-axial test device; a pressure transducer, coupling to the inflating fluid outlet, configured to measure pressure of the inflating fluids; and a computer, coupling to the LVDT and the pressure transducer, configured to gather and plot the measured pressure and vertical displacement throughout the bi-axial test.
 12. A method of numerically analyzing structure that contains elastomer comprising: conducting a bi-axial tension test of a specimen of the elastomer using a bi-axial test device of an elastomer material properties determination system, wherein the bi-axial tension test comprises loading, unloading and reloading phases; plotting measured pressures versus vertical displacements throughout the loading, unloading and reloading phases; determining a first set of material constants by a least square fitting technique using the measured pressures and displacements in the unloading phase; and determining a second set of material constants by least square fitting technique using the measured pressures and displacements in the reloading phase.
 13. The method of claim 12, wherein the specimen is a membrane made of the elastomer in uniform thickness.
 14. The method of claim 13, wherein the vertical displacement is at the center of the membrane.
 15. The method of claim 12, further comprising: numerically analyzing the structure using a constitutive equation that describes behaviors of the elastomer including the loading phase, the unloading phase and the reloading phase separately.
 16. The method of claim 15, wherein the constitutive equation includes Mullins effect of the elastomer.
 17. The method of claim 12, said conducting a bi-axial tension test further comprises placing the bi-axial test device in a temperature control environment.
 18. The method of claim 17, wherein the temperature control environment is a water bath or an environment controlled chamber. 